Effective mass suppression in a ferromagnetic two-dimensional electron liquid 
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We present numerical calculations of the electron effective mass in an interacting, ferromag- 
netic, two-dimensional electron system. We consider quantum interaction effects associated with 
the charge-density fluctuation induced many-body vertex corrections. Our theory, which is free of 
adjustable parameters, reveals that the effective mass is suppressed (relative to its band value) in 
the strong coupling limit, in good agreement with recent experimental results. 
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I. INTRODUCTION 

Two-dimensional electron systems (2DESs) realized 
at semiconductor interfaces are of continuing interest)^ 
from both basic physics and technological points of view. 
As a function of the interaction strength, which is char- 
acterized by the ratio of the Coulomb energy to Fermi 
energy, many novel correlated ground states have been 
predicted such as a paramagnetic liquid (rg < 26), fer- 
romagnetic liquid (26 < Ts < 35) and Wigner crystal 
(rs > 35)^. In the paramagnetic liquid phase, interac- 
tion typically leads to an enhancement of effective mass 
(to*) and spin susceptibility (x* cx g*m*), where g* is the 
Lande g*-factor. Effective mass is an important concept 
in Landau's Fermi liquid theory since it provides a direct 
measure of the many-body interactions in the electron 
system as characterized by increasing rs- 

The effective mass to* renormalized by interactions 
has been experimentally studied^'^'^i^i'^ for various para- 
magnetic 2DESs as a function of r^. In the highly in- 
teracting, dilute, paramagnetic regime (3 < Ts < 26), 
TO* is typically significantly enhanced compared to its 
band value, TOf,, and tends to increase with increasing 
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A question of particular interest 
is the dependence of to* on the 2D electrons' spin and 
valley degrees of freedom as these affect the exchange in- 
teraction. Recent measurements of to* for 2D electrons 
confined to AlAs quantum wells revealed that, when the 
2DES is fully valley- and spin-polarized, to* is suppressed 
down to values near or even slightly below TOa^ d^d^ . Note 
that in these experiments, < 22 so that the 2DES is 
in the paramagnetic regime, but a strong magnetic field 
is applied in order to fully spin-polarize the electrons. 
Here we present theoretical calculations indicating that 
the TO* suppression is caused by the absence (freezing 



out) of the spin fluctuations. The results of our to* cal- 
culations are indeed in semi-quantitative agreement with 
the measurements. 

Previous theoretical calculations of the effective mass 
are mostly performed within the framework of Landau's 
Fermi liquid theory whose key ingredient is the quasi- 
particle (QP) concept and its interactions. This entails 
the calculation of effective electron-electron interactions 
which enter the many-body formalism allowing the cal- 
culation of effective mass. A number of works consid- 
ered different variants of the leading order in the screened 
interaction for the self-energ y -^°did^d'^dSd9i20i2i,22 fj-Qj^ 

which density, spin-polarization, and temperature depen- 
dence of effective mass are obtained. In these calculations 
the on-shell approximatio n^^'^°i^^ yields a diverging ef- 
fective mass but the full solution of the Dyson equation 
yields only a mild enhancement. ^"'^^'^^ Almost all these 
works considered a paramagnetic 2DES as past experi- 
ments concentrated on the effective mass enhancement 
in partially spin polarized 2D systems with < 26. 



II. THEORY 

We consider a ferromagnetic 2DES as a model for a 
system of electronic carriers with band mass mi, in a 
semiconductor heterostructure with dielectric constant 
K. The bare electron-electron interaction is given by 
Vq = 27re^/(Kg). At zero temperature there is only 
one relevant parameter for the homogeneous, ferromag- 
netic 2DES, the usual Wigner-Seitz density parameter 
Ts = (Trna^)^^/^ in which as = /i^K/(TOf,e^) is the Bohr 
radius in the medium of interest. 

The QP self-energy with momentum k and frequency 
w in a fully polarized electron system can be written as 
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FIG. 1: (color online). Many-body effective mass as a func- 
tion of Ta for < < 22 for a ferromagnetic 2DES. 



(47rn2D)^^^ = 2/{rsaB), respectively, being the Fermi 
energy and wave vector; nF(fc) is the Fermi function. 
In Eq. Il]), e(q, is the dynamical screening function 
for which we use the form appropriate for a ferromag- 
netic 2DES derived from Kukkonen-Overhauser effective 
interaction.^- The many-body exchange and correlation 
(XC) effects are introduced through the local-field fac- 
tors (LFF) Gct.ct' (q, i^) (cr and a' are spin indices) take 
the Pauli-Coulomb hole around a charged particle into 
account. The dynamical screening function reads 
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l-G+(q,c^) xc(q,c^), (2) 



where is the LFF associated with charge-fluctuations. 
This expression is similar to the Kukkonen and Over- 
hauser interaction^^ where the spin-fluctuation term is 
dropped. A similar expression has also been reported in 
Refs. [23125^ . In Eq. (2) xc{<l,'^) represents the density- 
density response function, which in turn is determined 
by the local- field factor (q, uj) via the relation 
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in which x'j'('i',^) is the density response function of the 
spin-polarized electrons. The expression for the nonin- 
teracting density response function on the imaginary fre- 
quency axis is obtained for use in Eq. ^ as 




where we have defined flj = / 4ml — kp /m1 — Vi^/h^. 
It is evident that setting Gj'(q, (jj) = 0, we recover the 
standard random phase approximation (RPA). In what 
follows, we shall make the common approximation of ne- 
glecting the frequency dependence of . 

Quite generally, once the QP retarded self-energy is 
known, the QP excitation energy (5£Qp(k), which is the 
QP energy measured from the chemical potential /i^ of 
the interacting ferromagnetic 2DES, can be calculated by 
solving self-consistently the Dyson equation 



'54pW=^k+3?eSfet(k,^)| 



(5) 



Alternatively, the QP excitation energy can also be cal- 
culated from 



<54p(k) 



(6) 



This is called the on-shell approximation (OSA) and it 
is argued-- to be a better approach than solving the 
full Dyson equation since noninteracting Green func- 
tion is used in Eq. (1). Here 3fJeI]^gj (k, is defined as 



5ReSj,,(k,u;) 



The effective mass m*{k) is now calculated from 
1 1 d5£l.p{k) 



mUk) h^k 



dk 
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where for SE^p we have at our disposal the Dyson and 

OSA approaches. Evaluating m*(fc) at k = kp, one gets 
the QP effective mass at the Fermi contour. Clearly from 
Eqs. ([2]) and ([3]) LFF is the basic quantity for an evalu- 
ation of the QP properties. We have used the param- 
eterized forms of LFFs G~^{qX) and G~{qX) (and in 
particular G^{q) = G~^{q,( = 1) where ( is the spin 
polarization) of Moreno and Marinescui^ 



III. RESULTS AND DISCUSSION 

We now present our numerical results, which are based 
on the LFF G^{q) as input. In Fig. 1 we show our nu- 
merical results of the QP effective mass both in OSA 
and Dyson approximations. The QP effective mass sup- 
pression is substantially smaller in the Dyson equation 
calculation than in the OSA; the reason is that a sig- 
nificant cancellation occurs between the numerator and 
the denominator in the effective mass expression in the 
Dyson approach. To clarify the effect of charge-density 
fluctuation we have also shown the RPA results which do 
not take the strong many-body fluctuations into account. 
Note that the LFF takes into account multiple scattering 
events to infinite order as compared to the RPA where 
these effects are neglected. In the limit of small ( and 
rs 0, the effective mass can be analytically shown to 
be mymb = 1 -h (1 - C/2.0)r^ lnrs/(V27r) which our nu- 
merical calculations faithfully reproduce. 
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FIG. 2: (color online). Many-body effective mass as a func- 
tion of Ta for < r-j, < 22 for the ferromagnetic 2DES in 
comparison to experiments in Ref. [illl^. Different symbols 
denote different samples; triangles: A, squares: B, circles: C, 
and diamonds: D. 



In Fig. 2 we compare our effective mass calculations 
with the experimental results^'^^'^^. The measurements 
were made on 2DESs confined to modulation-doped AlAs 
quantum wells (QWs) of width 4.5, 11, 12, and 15 nm 
(samples A, B, C, and D). These samples were grown 
on GaAs substrates using molecular beam epitaxy. In 
bulk AlAs, electrons occupy three degenerate ellipsoidal 
conduction band valleys at the X-points of the Brillouin 
zone with longitudinal and transverse effective masses 
mi=1.05 and mf=0.205 (in units of the free electron 
mass). Thanks to the slightly larger lattice constant of 
AlAs compared to GaAs, the AlAs QW layer is under 
bi-axial compressive strain. Because of this compres- 
sion, the 2DES in the wider QW samples (B, C, and 
D) occupy two in-plane valleys with their major axes ly- 
ing in the plane^'*. In our measurements on these sam- 
ples, we applied uni-axial, in-plane strain to break the 
symmetry between these two valleys so that only one 
in-plane valley, with an anisotropic Fermi contour and 
band effective mass of mb — ^mirrit — 0.46 is occupied^^. 
In sample A, however, thanks to its very small QW 
width, the confinement energy of the out-of-plane val- 
ley is lower (because of its larger mass along the growth 
direction), so that the electrons occupy this valley and 
therefore have an isotropic Fermi contour and band ef- 
fective mass is mt = nit = 0.205^. The effective masses 
were deduced from the temperature dependence of the 
Shubnikov-de Haas oscillations, the details of which are 
given in Refs.— '^^i^^ . We emphasize that the data shown 
here (Fig. 2) were taken on single-valley 2DESs which 
were subjected to sufficiently large magnetic fields to fully 
spin polarize the electrons . 



It appears in Fig. 2 that the OSA accounts overall for 
the observed reduction of r7i| below the band value rea- 
sonably well. The agreement is particularly good for the 
wider samples (B, C, and D) which have > 7. The m* 
data for the narrowest sample (A), however, fall above 
the theoretical predictions. We do not know the reason 
for this discrepancy. However, we point out that, besides 
the difference in the shapes of the Fermi contour, there is 
another difference between sample A and the other three 
samples. Because of the very narrow width of sample A's 
quantum well and the prevalence of interface roughness 
scattering^^, the mobility of the electrons in this sample 
is much lower (about a factor of 6) than in other sam- 
ples for comparable rg. It is possible that the higher 
disorder in sample A is responsible for m* being larger; 
this conjecture is indeed consistent with the results of 
calculation^ which predict a larger m* for more disor- 
dered samples. 

From Figs. 1 and 2 we draw two main conclusions, (i) 
The RPA and present results are rather similar in the 
weak coupling limit (ts < 1). (ii) In the strong coupling 
regime (r^ > 3), however, our theoretical calculations 
which incorporate the proper many-body effects exhibit a 
mass suppression, similar to the experimental data, while 
the RPA results show a mass enhancement and are far 
from the experimental data. We emphasize that the effec- 
tive mass at the Fermi contour is significantly suppressed 
in the fully polarized case because of the absence of spin- 
fluctuation contribution. This suppression suggests that 
the anti-symmetric Landau parameter Fi < and thus 
higher angular momentum Landau parameters may be 
negligible in a fully spin-polarized 2DES. 
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FIG. 3: (color online). Many-body on-shell effective mass as 
a function oi k/kp at rs = 5 for 2DES with the combined 
effect of charge fluctuations in comparison to paramagnetic 
2DES. 

To gain further insight to the density dependence of 
m* , we have calculated the on-shell effective mass as a 
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FIG. 4: (color online). Renormalization constant as a 
function of for < < 22 for a ferromagnetic 2DES. 



function of particle momentum fc using Eq. ([7]) evaluated 
at = 'Ck/?! and rs = 5. More specifically, we use 



(8) 



for a ferromagnetic case. The results for both paramag- 
netic and ferromagnetic cases are shown in Fig. 3. m*(fc) 
for a paramagnetic 2DES by using G^{q,C = 0) and 
G~(g, C ~ 0) has a sharp peak around k kp where 
fcp = 2/(rsaB) and a resonance like divergent behavior 
around k sa 2kp. The peak around fcp is associated with 
spin fluctuations and the divergent behavior around 2fcF 
is related to density fluctuations.— In particular, the 
latter divergence has been extensively studied by Zhang 
et al^ within the RPA. It is related to the dispersion in- 
stability and coincides with the plasmon emission, m^(fc) 
for the ferromagnetic 2DES, on the other hand, clearly 
shows the disappearance of the peak associated with spin 
fluctuations. Thus, TO*(fc) is very weakly momentum de- 
pendent for k < kp, since there is a substantial cance- 
lation between the residue and the exchange plus line 
self-energy contribution in this regime which make the 
real part of the retarded self-energy approximately lin- 
ear with respect to fcjii The divergence associated with 
charge fluctuations is still present, showing a negative 
peak around k = 2kp. m|(fc) calculated within the RPA 
reproduces quantitatively the divergent behavior associ- 
ated with charge fluctuations but shows some structure 
for k < kp, therefore failing to account for the absence of 
spin fluctuations. Our density-dependent effective mass 
results (Figs. 1 and 2) are consistent with m*{k) calcula- 
tions which we have checked for a range of values. 



We have also calculated the renormalization factor 
Z|(rs) which is equal to the discontinuity in the mo- 
mentum distribution at kp and defined by Z^^ = 

1 — h^^ duj'3ieT,l^^{k,uj) . The effect of charge- 

k—kp 

fluctuations is to make the Z| values larger at large 
compared to the case when they are not included as 
shown in Fig. 4. This means that charge-density fluc- 
tuations tend to stabilize the system, whereas the RPA 
works in the opposite direction.— In the present case in- 
cluding the LFF helps preserve the Fermi liquid picture 
in the low density regime. 

We have performed our numerical calculations for 
strictly 2DES. As indicated above, experimental samples 
have a finite thickness in the range of 5—15 nm. Our theo- 
retical model may be extended to include the finite quan- 
tum well width effects in the following manner. Choos- 
ing, say, an infinite square well model with width L will 
modify the bare Coulomb interaction, Vq — > VqF{qL), 
in which F{x) is a form factor.— For consistency, one 
should also calculate the local-field factor G^{q) using 
the same model for the finite width effects. This would 
provide a better comparison with experiments. In our 
case, the local-field factor we use was constructed22. by 
the quantum Monte Carlo data for a strictly 2DES and 
it is not straightforward to incorporate the finite with ef- 
fects within such an approach. Previous calculations^^ of 
the effective mass for a paramagnetic 2DES suggest that 
the effect of a finite thickness is to suppress m* . There- 
fore, we surmise that a similar qualitative effect would 
occur for the ferromagnetic 2DES. On the other hand, 
the finite temperature and disorder effects have a ten- 
dency to enhance the effective mas o^°i^^ which may lead 
to a cancellation. These issues require a more systematic 
study. 



IV. SUMMARY 

In conclusion, our theoretical calculations incorporat- 
ing the proper Pauli-Coulomb hole and multi-scattering 
processes show that in an interacting, fully spin-polarized 
2DES the absence of spin fluctuations reduces the effec- 
tive mass below its band value, in agreement with exper- 
imental data. Our results also demonstrate the inade- 
quacy of RPA to account for the observed effective mass 
suppression. 
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